Dynamical Mean Field Theory of Correlated Gap Formation in Pu Monochalcogenides 
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The correlated Kondo insulator state of the plutonium monochalcogenides is investigated using 
the dynamical mean field theory (DMFT) and the local density approximation +[/ (LDA+[/). The 
DMFT-dynamical fluctuations lead to a correlated insulator state at elevated temperature, in sharp 
contrast to the static LDA+[/ approach that fails to reproduce both the insulating behavior and 
anomalous lattice constant. The DMFT conversely predicts the experimentally observed anomalous 
increase of the gap with pressure and explains the lattice constant very well. 

PACS numbers: 71.20.Nr, 71.27.+a, 75.30.Mb 



Correlated /-electron materials pose a major challenge 
to electronic structure theories. In particular, plutonium- 
based materials have received considerable attention re- 
cently because of the puzzling complexity of their cor- 
related 5/ electron behavior, for which so far no unified 
theoretical treatment could be established. The com- 
plex nature of Pu materials is reflected in the variety 
of widely different ground states that may emerge out 
of the competition of electronic and magnetic interac- 
tions involving the Pu 5/ electrons. For instance, the 
non-magnetic ground state of the open 5/ shell of J-Pu 
is a long standing controversial issue, and the excep- 
tional occurrence of superconductivity at a very high crit- 
ical temperature of 18.5 K discovered in PuCoGas and 
its relation to the 5/ electronic structure is an unsolved 
mystery^^ A correlated insulating state that is very simi- 
lar to that of Kondo insulators has been observed for the 
Pu monochalcogenides.^^^ 

Energy gaps in the excitation spectrum of corre- 
lated /-electron materials, such as SmBe and Ce3Bi4Pt3, 
have been intensively investigated in the last years.^ 
These Kondo insulator materials typically exhibit heavy- 
fermion characteristics in conjunction with narrow gap 
semiconducting behavior.^ The quasi-particle gap often 
shows an unusual dependence with pressure or tem- 
perature, which indicates the existence of a narrow /- 
derived feature near the top of the valence band. The 
Pu monochalcogenides show many anomalous proper- 
ties that are paramount to those of Kondo insulators. 
In particular, they are narrow gap semiconductors with 
an unusual temperature and pressure dependence of 
the gap,-i^ yet the specific heat is surprisingly high, 
7 = 30 mJ/molK^.^ Although electronic structure in- 
vestigations of the Pu monochalcogenides have been 
performed, ^^^^^^^^ these could not provide a consistent 
explanation of the anomalous properties. 

Here we show that the anomalous properties of the 
Pu monochalcogenides can be well explained by dy- 
namical mean field theory (LDA+DMFT) calculations. 
In particular, the non-magnetic ground state, the for- 
mation of a correlated gap at moderate temperatures, 
the anomalous increase of the gap with pressure and 
the unusual equilibrium lattice parameter are precisely 
given by LDA+DMFT theory. Remarkably, although 



the DMFT is conventionally considered to be ideally 
suited for correlated materials with a sharp quasi-particle 
resonance occurring at the Fermi energy, ^^^^^ we report 
here the formation of a Kondo insulator gap within the 
LDA+DMFT. We compare furthermore the results of 
DMFT and LDA+[/ calculations. Although the LDA+/7 
is considered to be the static limit of the DMFT, we find 
that the LDA+/7 approach /ai/5 to describe the correlated 
gap formation as well as the anomalous lattice constant, 
emphasizing the importance of dynamical fluctuations for 
the quasi-particle gap formation. 

The physical properties of the Pu monochalco- 
genides (PuS, PuSe, and PuTe) are experimentally 
well documented through transport, susceptibility, 
specific heat, and photoelectron spectroscopy (PES) 
measurement s.^i^i^ii^ii^ The lattice constants of the Pu 
monochalcogenides are anomalous, as they are too small 
for both divalent (5/^) and trivalent (5/^) Pu ions, which 
signals that these actinide chalcogenides are not plain 
ionic salts. Magnetic order, furthermore, is absent 1^ The 
temperature dependence of the resistivity signals a com- 
plex semiconducting behavior with temperature depen- 
dent gaps.^ Under pressure the room-temperature gap in- 
creases completely contrary to conventional knowledge.^- 
The specific heat coefficient is quite high and difficult to 
reconcile with the semiconducting behavior. Photoemis- 
sion experiments performed on PuSe and PuTe reveal a 
sharp /-related peak near the Fermi edge^i^ii^ 

The peculiar phenomena of the Pu monochalcogenides 
have been attributed to intermediate valence behavior by 
Wachter et al.^ In this explanation the gapping near the 
Fermi energy is modeled through adopting a special form 
of the density of states (DOS). LDA calculations predict 
a low DOS at the Fermi level, but fail to predict the 
insulator state and also predict a too small equilibrium 
lattice parameter (Ref. 10 and references therein). Multi- 
plet calculations within the LDA-Hubbard I approxima- 
tion (LDA-HIA) were recently used— to successfully 
explain the peak structure present in the experimental 
PES spectra. Also, an LDA+/7 and LDA+DMFT inves- 
tigation of actinide monochalcogenides was recently per- 
formed by Pourovskii et al}^ These calculations predict 
an almost completely filled 5/5/2 subshell (i.e., divalent 
state) for the Pu monochalcogenides and do not repro- 
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duce an energy gap. Consequently, there is presently no 
theory available to suitably explain the mechanism driv- 
ing the correlated gap formation of these Kondo insula- 
tors. We focus here on the development of this insulating 
state, as well as the unusual increase of the gap under 
pressure. 

We performed electronic structure calculations using 
the accurate, relativistic full-potential, linearized aug- 
mented plane wave (FLAPW) method, in which the spin- 
orbit interaction is implemented through a second vari- 
ational step. For the LDA+/7 calculations we used the 
around mean field form of the double counting term. We 
adopt here for the screened Coulomb and exchange pa- 
rameters U = 3 — 4 eV and J = 0.55 eV, which are in the 
range of commonly accepted values for Pu compounds. 
The two-electron integrals of the Coulomb interaction of 
the / electrons are expressed in terms of effective Slater 
integrals F^, (0 < < 6), where Fq = [/, and, in a man- 
ner similar to that of Ref. 2, we rescale the parameters. 

The ideas of the LDA+DMFT method have been re- 
viewed previously ] ^^i^^i-*^^ We compute the local Green 
function employing the Kohn-Sham states obtained from 
the relativistic LDA+/7 calculation in a manner simi- 
lar to Ref. 19. The local Green function is computed on 
Matsubara frequencies uj^ through a Brillouin zone (BZ) 
integration of the lattice Green function, Gbb'ik^ioOn) = 
{bk\[{iujn + M — Sbk)^ — AT^{iuJn)]~^\b'k), where /i is the 
chemical potential, bk the lattice states, and AT^{iuJn) the 
local self-energy. The latter is corrected for the static 
Hartree-Fock self-energy part that is already contained 
in the LDA+/7 calculation, i.e., AT^{iUn) = T^{iuJn) — 
Shf- The spin-polarized T-matrix fiuctuation-exchange 
(FLEX) impurity solver^^^^^ is used here for generating 
the self-energy. We note that the FLEX impurity solver 
is expected to be applicable to moderately correlated 
materials, but it cannot sufficiently resolve the atomic 
multiplet structure. This implies that with the present 
DMFT-FLEX formulation the PES spectrum, containing 
atomic multiplets,^^'^^ cannot be fully captured. How- 
ever, DMFT-FLEX is precisely suitable for our objective, 
namely, the investigation of the correlated gap formation 
at elevated temperatures. 

In our calculations we used 231 /c-points in the irre- 
ducible BZ for the self-consistent LDA+/7 calculation 
and the DMFT self-consistency loop. In the later we 
take into account the / orbit als and their spin, projected 
from all energy bands, resulting in a 14 x 14 matrix for 
the local Green function and self-energy. Overall self- 
consistency has been achieved through iterative feeding 
of the density matrix of the local Green function in the 
next LDA+/7 loop, and back feeding of the new solutions 
in the DMFT loop. The /c-dependent and /c-integrated 
DOS was computed on 1469 /c-points, using a Pade ap- 
proximation to G{k^iuj). In the FLAPW calculation we 
included the Pu 5/, 6(i, 75, and 7p orbitals as band states 
and the Pu 5(i, 65, and 6p states as semi-core states, 
whereas for the chalcogenides we treated the ns, np, and 
nd states as band states and the {n — l)d states as semi- 
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FIG. 1: (Color online) Top: the energy bands of PuTe com- 
puted with the LDA+[/ method at the experimental lattice 
constant. Bottom: LDA+DMFT quasi-particle band struc- 
ture of PuTe at T = 600 K. The colors indicate the magnitude 
of the /c-dependent spectral function, — 7r~^Im G{k,E). 

core states, as well as the {n — l)p states for PuTe. In 
the DMFT calculation, we furthermore used 8192 Mat- 
subara frequency points and applied an analytic asymp- 
totic form of the Green function to compensate for the 
remaining Matsubara frequencies when taking the sum 
over the frequencies. Although we use a large number 
of Matsubara frequencies, we can only compute a finite 
number of frequency points and therefore the calcula- 
tions are restricted to moderately high temperatures in 
practice (300 K and above). This corresponds precisely 
to the temperature range in which the room temperature 
gap of the Pu monochalcogenides has been measured. 

In Fig. [U we show the energy band structure of PuTe 
calculated with the LDA+/7 method at the experimen- 
tal lattice constant. As was noted before, the around 
mean field LDA+/7 approach (as well as localized limit 
LDA+[/) increases the spin-orbit interaction, favors j-j 
coupling and hence drives many Pu systems to a non- 
magnetic ground state. ^^^^ For the Pu monochalcogenides 
this agrees nicely with the experimentally detected ab- 
sence of magnetic ordering.^ However, as Fig. [1] reveals 
the LDA+/7 method predicts the Pu monochalcogenides 
to be semi- metallic. Although the DOS at the Fermi level 
is not high, the obtained metallicity is in contradiction 
with experimental knowledge. This finding gives a first 
indication that the energy gap of the Pu monochalco- 
genides is special. For other correlated actinide com- 
pounds, such as UO2 and PUO2, the standard LDA ap- 
proach fails to predict an insulator state, but the LDA+/7 
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FIG. 2: (Color online) The /c- integrated DOS calculated with 
the LDA+DMFT method for (a) PuTe and (b) PuSe. The 
experimental photoemission spectraF^?^ are given by the solid 
and dashed black lines. 



approach, which is designed to give the static self-energy 
correction to localized states, enforces the formation of 
a lower and an upper Hubbard band. However, this 
typical Mott-Hubbard scenario does not apply to the 
Pu monochalcogenides. In Fig. [T] we also show the k- 
dependent spectral density of PuTe computed with the 
LDA+DMFT. Notably, the LDA+DMFT predicts the 
formation of a narrow gap in the quasi-particle spectrum. 
A detailed comparison of the LDA+[/ and LDA+DMFT 
bands in Fig. [1] provides insight in the influence of the 
dynamical fluctuations in the self-energy that form the 
difference between the static LDA+/7 and DMFT. The 
dispersive energy bands that occur at binding energies of 
2 to 6 eV are dominantly due to the telluride bp states 
with a small admixture of Pu 5/ states; these bands are 
obviously unaffected. The narrow, / derived band below 
the chemical potential changes more: the LDA+/7 calcu- 
lations shift this hybridized / state to a lower energy, as 
compared to LDA calculations. Consequently, the 5/ oc- 
cupation is increased in the LDA+/7. In the DMFT, the 
fluctuations within the / shell are dynamically screened 
by fast conduction electrons, which leads to an upward 
shift of the narrow state. The hybridization of the / 
states and the conduction band near the chemical poten- 
tial is increased, which causes the formation of a small 
energy gap. We note that the DMFT calculation, start- 
ing from the non- magnetic LDA+/7 result, does not lead 
to a magnetic ground state. Also, the DMFT reduces the 
/ occupation with regard to that of the static LDA+/7. 

In Fig. [2] we show the /c-integrated DOS obtained with 
the LDA+DMFT method for PuSe and PuTe. As most of 
the recent experiments have been performed on PuSe and 
PuTe^ii^ii^ we concentrate on these two chalcogenides. 
The experimental PES spectra^^^^^ are also shown in Fig. 
[2] for comparison. These experimental data contain the 
"three-peak" structure at energies between and —1 eV, 
followed by a broad peak at — 2 eV and a second, smaller 
hump at —4 eV. Considering first PuSe, we recognize 
that the second broad peak at — 4 eV is predicted to be 
due to the Se p-states. This finding is consistent with 
He I and He II PES^^ which attributes the three-peak 
structure and the peak at — 2 eV to /-related excitations, 
but not the hump at —4 eV. Our DMFT calculations 
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FIG. 3: (Color online) The relative total energy of PuTe ver- 
sus lattice constant, calculated using the LDA, LDA+t/ and 
LDA+DMFT approaches. The experimental value is indi- 
cated by the vertical dashed line. Open circles and crosses 
indicate metallic and insulating ground states, respectively. 

predict a narrow /-related peak just below the Fermi 
edge for PuSe and PuTe, which corresponds to a steep 
peak in the experimental spectra. The /-related peaks at 
—2, —1 and —0.5 eV have previously been attributed to 
multiplet structuresi2iiii and these are not captured by 
our DMFT-FLEX calculations. Interestingly, the steep 
peak closest to the Fermi edge has been attributed to a 
multiplet, too^iSiiii As a salient difference, our DMFT- 
FLEX calculations predict this peak instead to be due 
to a narrow, hybridized / band (see Fig. [1]). Wachter^^. 
used the magnetic susceptibility to argue that the peak 
near the Fermi level should be due to a narrow 5/ band, 
and not a multiplet peak. We propose that temperature- 
dependent PES measurements could be used to discrimi- 
nate between these two explanations, as multiplet peaks 
are not expected to show a temperature dependence. 

The integrated DOS of PuSe and PuTe were also com- 
puted by Pourovskii et a/.,^^ however, we obtain here a 
DMFT solution that is quite different. Pourovskii et ali^ 
computed a chalcogenide DOS that is shifted towards 
the Fermi edge by a considerable 2 eV as compared to 
our results, and hence, they attribute both the exper- 
imental peaks at —2 and —4 eV to the chalcogenide p 
DOS. In addition, they obtain a Pu configuration near 
the closed shell divalent state, with an / occupation of 
5.8, both in their LDA+/7 and non-selfconsistent (single- 
shot) DMFT calculations. Our LDA+/7 implementation 
conversely gives an intermediate valence 5/ occupation 
of about 5.5, which, in the self-consistent LDA+DMFT 
loop is reduced to 5.3, as the / state is depleted when it 
energetically comes closer to fi. Thus, our DMFT calcu- 
lation gives still an intermediate valence occupation, but 
quite close to the trivalent state, whereas the previous 
study^^ predicted it to be almost a 5/^ configuration. 
Other recent DMFT studies^^^^ for ^-Pu also computed 
Pu configurations relatively close to 5/^. 

The small energy gap predicted by the DMFT could 
be a fortuitous coincidence. We show in the following 
that this is not the case. First, in Fig. [3] we plot the to- 
tal energy of PuTe versus lattice parameter a, computed 
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FIG. 4: (Color online) (a) Pressure dependence of the energy 
gap of PuTe. The blue curve shows the experimental gap^ and 
the red curve depicts the LDA+DMFT result, (b) The quasi- 
particle bands around the chemical potential for two lattice 
constants, corresponding, respectively, to ambient pressure 
and -5.38 GPa. 

with the LDA, LDA+/7, and LDA+DMFT^"^ approaches. 
The LDA predicts a metalHc ground state and a much 
too small equilibrium lattice parameter (cf. Ref.ilO|). The 
LDA+/7 also predicts a metallic ground state, but with 
a too large aS^ Only the LDA+DMFT predicts the in- 
sulator state and the experimental lattice constant. 

Lastly, we discuss the pressure dependence of the cor- 
related energy gap, which is a further stringent test of 
the DMFT solution. We find that pressure causes a 
stronger mixing of Pu 5/ and conduction states near 
the Fermi level, leading to a widening of the gap. The 
experimental^ and computed gaps of PuTe as a function 
of pressure are presented in Fig. 2^. The DMFT gap val- 
ues were computed from the maximum intensity points 
of the highest occupied and the lowest excited state. The 



experiment reported that the room-temperature activa- 
tion gap first increases with pressure up to 5 GPa and 
then slightly decreases.^ Our DMFT calculations com- 
pletely reproduce and explain this. The quasi-particle 
bands, shown in Fig. [4)3, have the highest occupied state 
at the r point. The lowest excited states are along the 
r — X and U — X directions for low pressures, but at 
the F point for pressures above 5.38 GPa. Thus, our cal- 
culations predict that under pressure a subtle (indirect- 
direct) transition in the /c-dependent lowest excited states 
occurs, which gives rise to a turn in the slope of the 
gap versus pressure curve. The gap increase is counter- 
intuitively, however, it was shown that either a gap in- 
crease or decrease can occur for Kondo insulators. 

Previously there existed no explanation of the peculiar 
insulator state of the Pu monochalcogenides; electronic 
structure approaches, as the LDA and LDA+[/ that are 
based on static approximations to the self-energy, fail to 
reproduce the insulator state. We have identified here 
that dynamical fluctuations in the DMFT are the source 
that drives the formation of the correlated insulator state. 
Our LDA+DMFT calculations accurately explain the 
lattice constant as well as the anomalous increase of the 
room- temperature gap under pressure. Also, we find that 
the Pu 5/ occupation, which is given to be intermediate 
between 5 and 6 by the LDA+/7 approach, is reduced in 
our LDA+DMFT calculations, and consequently, comes 
closer to the trivalent configuration. Our LDA+DMFT 
calculations predict furthermore the presence of a narrow 
Pu 5/-derived band close to the Fermi edge. We propose 
that the existence of the narrow-band state could be ver- 
ified through temperature-dependent PES experiments. 
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